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ABSTRACT 


The  structural  dynamic  response  of  turbomachinery  components 
to  aer 6-thermodynamic  distortion— induced  excitations  is  an  item 
of  major  concern  in  the  design  of  advanced  gas  turbine  engines. 
The  rotor  speeds  at  which  these  resonant  forced  responses  occur 
can  be  predicted  with  Campbell  diagrams.  However,  due  to  the 
inadequacies  of  the  existing  time-variant  aerodynamic  models,  no 
accurate  prediction  can  currently  be  made  for  the  amplitude  of 
the  resulting  stresses. 

.  ,  0 

The  overall  objective  of  this  research  program  is  to  quanti¬ 
tatively  investigate  the  fundamental  phenomena  relevant  to  aero- 
thermodynamic  distortion  —induced  structural  dynamic  blade 
response  in  multi-stage  gas  turbine  fans  and  compressors.  Unique 
unsteady  aerodynamic  data  will  be  obtained  to  validate  and  indi¬ 
cate  necessary  refinements  to  state-of-the-art  analyses  and  to 
direct  the  modeling  of  new  analyses.  Also,  for  the  first  time,  a 
first  principles  capability  to  predict  the  vibrational  response 
amplitude  of  blading  due  to  aerodynamic  excitations  will  be 
developed.  - 

\ 

This  report  summarizes  the  "progress  and  results  obtained 
during  the  first  year  of  this  program^  Thesej? include:  the 
dynamic  instrumentation  of  the  first-stage  vane  row  of  a  three- 


stage  research  compressor;  the  design,  specification,  and  initia¬ 
tion  of  the  development  of  the  dynamic  data  acquisition  and 
analysis  system;  the  initiation  of  the  development  of  the  neces¬ 
sary  data  analysis  and  calibration  techniques;  and  the  theoreti¬ 
cal  development  of  a  unique  coupled  mode  structural  dynamic  blade 
response  analysis  based  on  an  energy  balance  technique. 
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1 .  INTRODUCTION 


The  structural  dynamic  response  of  fan,  compressor,  and  tur¬ 
bine  blading  to  aero-thermodynamic  distortion  induced  excitations 
is  an  item  of  rapidly  increasing  concern  to  designers  and 
manufacturers  of  gas  turbine  engines  for  advanced  technology 
applications.  Destructive  aerodynamic  forced  responses  of  gas 
turbine  engine  blading  have  been  generated  by  a  wide  variety  of 
aero-thermodynamic  distortion  sources.  These  include:  blade 
wakes;  multi-staging  interaction  effects;  large  angle  of  attack 
or  yaw;  engine  exhaust  recirculation;  cross-flow  at  the  inlet; 
pressure  variations  on  the  engine  due  to  external  aerodynamics; 
and  armament  fixing.  The  following  physical  phenomena  have  been 
identified  as  being  significant  with  regard  to  aero-thermodynamic 
distortion  induced  structural  dynamic  responses  of  fan  and 
compressor  blading. 

*  Resonance 

*  Multi-stage  interactions 

*  Stall 

*  Inlet  gusts 

*  Time- varying  inlet  flows 


Flow  separation 


* 

*  Potential  effects 

*  Turbulence 


The  first  principles  prediction  of  the  structural  dynamic 
response  associated  with  all  of  the  above  physical  phenomena  is 
identical  and  involves  the  following  elements.  spatially 
periodic  variations  in  pressure ,  velocity,  temperature,  and  flow 
direction  of  the  exit  flow  field  of  an  upstream  element  appear  as 
temporally  varying  in  a  coordinate  system  fixed  to  the  downstream 
row.  As  a  result,  individual  airfoils  are  subject  to  a  time- 
variant  aerodynamic  forcing  function  which  can  induce  high  level 
vibratory  stresses. 

The  analysis  of  the  aerodynamically  forced  response  vibra¬ 
tory  behavior  of  a  blade  or  vane  row  requires  a  definition  of  the 
unsteady  forcing  function  in  terms  of  its  harmonics.  The  time- 
variant  aerodynamic  response  of  the  airfoil  to  each  harmonic  of 
this  forcing  function  is  then  assumed  to  be  comprised  of  two 
parts.  One  is  due  to  the  disturbance  being  swept  past  the  non¬ 
responding  fixed  airfoils.  The  second  arises  when  the  airfoils 
respond  to  this  disturbance.  Mathematically  these  effects  are 
modeled  by  two  distinct  analyses.  A  linearized  small  perturba¬ 
tion  gust  analysis  is  used  to  predict  the  time-variant  aerodynam¬ 
ics  of  the  fixed  non-responding  airfoils  to  each  harmonic  of  the 
disturbance.  A  self-induced  unsteady  aerodynamic  analysis 
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wherein  the  airfoils  are  assumed  to  be  harmonically  oscillating 
is  then  used  to  predict  the  additional  aerodynamic  effects  due  to 
the  airfoil  response.  Superposition  of  these  two  effects  can  be 
performed  only  with  knowledge  of  the  modal  pattern  and  amplitude 
of  response  of  the  blading  because  the  magnitude  of  the  pressure 
field  resulting  from  the  airfoil  motion  is  dependent  upon  the 
amplitude  of  this  motion.  Thus  a  solution  with  key  elements  con¬ 
sisting  of  the  gust  analysis,  a  self-induced  unsteady  aerodynamic 
analysis,  and  an  airfoil  structural  analysis  is  necessary  to 
predict  the  total  response  of  an  airfoil  subjected  to  an  upstream 
generated  spatially  periodic  disturbance. 

The  unsteady  small  perturbation  gust  and  self-induced 
unsteady  aerodynamic  analyses  are  two-dimensional  and ,  as  such, 
are  coupled  to  the  airfoil  structural  analysis  by  means  of  a 
strip  theory  approximation.  Thus,  the  airfoil  is  considered  to 
consist  of  a  series  of  individual  and  independent  two-dimensional 
aerodynamic  regions.  The  time-variant  aerodynamic  analyses  are 
then  applied  to  each  such  individual  region,  with  the  charac¬ 
teristic  parameters  including  the  Mach  number,  reduced  frequency, 
stagger  angle,  and  solidity,  taken  as  the  average  value  at  the 
inlet  boundary.  It  should  be  noted  that  there  is  no  coupling 
between  adjacent  aerodynamic  regions.  Hence,  aerodynamic  forced 
vibrations  involving  spanwise  variations  in  unsteady  aerodynamics 
cannot  be  treated,  i.e.,  if  spanwise  variatins  exist  in  the  aero¬ 
dynamic  forcing  f*  notion  ot  caused  by  simple  inlet  spanwise 
variatins  in  the  Mach  ..umber,  reduced  frequency,  stagger  angle, 


or  solidity,  such  a  strip  theory  design  system  is  of  no  value. 
Some  of  the  previously  noted  forced  response  sources  may  fit  into 
this  category,  as  for  example,  the  case  of  rotor  tip  vortices 
generating  a  forced  response  in  a  downstream  blade  or  vane  row. 


II.  PROGRAM  OBJECTIVES 


Overall  Objectives 

The  overall  objective  of  this  research  program  is  to  quanti¬ 
tatively  investigate  the  fundamental  phenomena  relevant  to  dis¬ 
tortion  generated  aero-thermodynamic  induced  structural  dynamic 
effects  in  gas  turbine  engine  blade  rows.  Unique  unsteady  aero¬ 
dynamic  data  will  be  obtained  on  stationary  and  rotating  blade 
rows  to  validate  and  indicate  necessary  refinements  to  current 
state-of-the-art  analyses  and  to  direct  the  modeling  and  develop¬ 
ment  of  new  analyses.  Also,  for  the  first  time,  a  first  princi¬ 
ples  capability  to  predict  the  vibrational  response  amplitude  of 
blading  due  to  aerodynamic  excitations  based  on  am  energy  balance 
technique  will  be  developed. 

From  first  principles  considerations,  the  relevant  fundamen¬ 
tal  physical  phenomena  are  identical  for  the  various  sources  of 
aero-thermodynamic  distortion.  Hence,  to  accomplish  this  overall 
objective  in  a  timely  and  efficient  manner,  while  obtaining 
results  of  direct  interest  and  significance  to  the  gas  turbine 
engine  community,  this  proposed  research  program  is  concerned 
with  the  time-variant  aerodynamics  and  structural  dynamic 
response  of  multi-stage  stationary  and  rotating  blade  rows,  with 
the  primary  source  of  excitation  initially  being  the  wakes  from 
upstream  blade  elements.  The  overall  objective  includes  the  fol¬ 
lowing. 


The  experimental  determination  of  the 
fundamental  time-variant  gust  aerodynamics 
associated  with  variations  in  incidence 
angle  (loading),  reduced  frequency, 
solidity,  and  multi-stage  effects  on 
both  stationary  and  rotating  blade  rows 
as  well  as  the  investigation  of  the 
validity  of: 

*  the  two-dimensional  linearity  and 
super positioning  assumptions; 

*  the  small  perturbation  modeling  concept; 

The  development  of  a  first  principles 
state-of-the-art  capability  to  predict  the 
vibrational  response  amplitude  of  blading 
based  on  an  energy  balance  technique,  with 
key  elements  consisting  of  the  following 
coupled  existing  analyses: 

*  a  flat  plate  cascade  transverse  gust 
analysis; 

*  a  flat  plate  cascade  self- induced 
unsteady  aerodynamics  analysis; 

*  an  airfoil  structural  analysis. 


Thus ,  this  progxam  is  directed  at  providing  fundamental 
time-variant  aerodynamic  data  which  will  not  only  validate  and 
indicate  refinements  to  the  current  state-of-the-art  two¬ 


dimensional  gust  analyses  but  will  also  address  the  validity  of 
the  most  basic  assumptions  inherent  in  these  analyses  and  in  the 
structure  of  forced  response  design  systems.  In  addition  a 
state-of-the-art  first  prinicples  predictive  aerodynamically 
forced  response  analysis  based  on  an  energy  balance  technique 
will  be  developed. 


Specific  Objectives 

The  following  were  the  specific  objectives  for  the  first 
year  of  this  research  program. 

*  Dynamic  instrumentation  of  the  first  stage 
vane  row  of  the  three-stage  research 
compressor; 

*  Design,  specification,  and  initiation  of  the 
development  of  the  dynamic  data  acquisition 
system; 

*  Initiation  of  the  development  of  the 
necessary  data  analysis  and  calibration 
techniques; 
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*  Theoretical  development  of  a  coupled  mode 
structural  dynamic  blade  response  analysis 
based  on  an  energy  balance  technique. 
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III.  TECHNICAL  APPROACH 


The  technical  approach  to  achieve  the  overall  program  objec¬ 
tives  requires  that  extensive,  high-quality,  detailed  aerodynamic 
data  be  acquired  from  benchmark  experiments  which  model  the  fun¬ 
damental  flow  physics  of  aero-thermodynamic  induced  structural 
dynamic  effects  in  gas  turbine  engines.  These  data  must  be 
analyzed  and  correlated  with  appropriate  state-of-the-art  ana¬ 
lyses  to  ascertain  their  range  of  validity  and/or  to  indicate 
necessary  refinements.  In  addition,  advanced  mathematical  models 
and  techniques  for  the  prediction  of  these  phenomena  will  be 
developed  and  experimentally  verified. 

The  approach  to  achieving  these  experimental  objectives  is 
to  measure  the  time-variant  pressure  distribution  in  controlled 
benchmark  experiments  which  model  the  fundamental  flow  physics 
and  thereby  identify  and  quantify  the  key  unsteady  aerodynamic 
parameters  relevant  to  aero-thermodynamic  distortion  induced 
response  of  blading.  In  particular,  this  research  program  is 
aimed  at  measuring  the  fundamental  unsteady  gust  aerodynamics 
associated  with  variations  in  incidence  angle  (loading),  reduced 
frequency,  solidity,  multi-stage  effects,  as  well  as  investigat¬ 
ing  the  validity  of  the  two-dimensional  linearity  and  superposi- 
tioning  assumptions  and  the  applicability  of  the  small  perturba¬ 
tion  modeling  concept.  This  requires  obtaining  airfoil  surface 
chordwise  and  spanwise  unsteady  pressure  data  on  both  stationary 


and  rotating  blade  cows  as  well  as  detailed  inlet  and  exit  plane 
flow  surveys  in  a  multi-stage  axial  flow  compressor.  These  data 
will  be  correlated  with  state-of-the-art  unsteady  cascade 
transverse  gust  analyses  and  will  be  used  to  determine  their 
range  of  validity  and  to  direct  the  development  of  an  advanced 
first  principles  forced  vibration  aerodynamic  analysis  which  the 
designer  can  use  with  confidence.  In  addition,  a  unique  first 
principles  predictive  aerodynamically  forced  response  analysis 
will  be  developed  based  on  an  energy  balance  technique. 


In  the  first  stage  of  a  multi-stage  axial  flow  compressor 
the  wakes  from  the  upstream  rotor  blades  are  the  source  of  the 
aerodynamically  induced  fluctuating  surface  pressures  on  the 
downstream  stator  vanes,  i.e.,  the  rotor  blade  wakes  define  the 
forcing  function  to  the  downstream  stator  vane  row.  In  the  aft 
stages  of  the  compressor  the  wakes  from  all  of  the  upstream  blad¬ 
ing,  both  rotors  and  stators,  are  the  source  of  the  aerodynami¬ 
cally  induced  fluctuating  surface  pressure  distributions  on  the 
downstream  stator  vane  and  rotor  blade  rows.  Hence,  it  is  neces¬ 
sary  to  expez imentally  model  all  of  the  significant  features 
which  define  this  forcing  function.  These  include  the  variation 
of  incidence,  the  wave  form,  the  velocity  (pressure)  variation, 
the  various  blade  row  interactions  in  a  multi-stage  machine,  and 
the  reduced  frequency  (k  -  wC/  2Vav t M , )  which  defines  the  order 
of  magnitude  of  the  unsteadiness.  (Forced  response  problems  in 


turbomachines  typically  have  first  harmonic  reduced  frequency 
values  varying  from  0.2  to  8.0.)  The  above  noted  features  can 
all  be  simulated  in  the  Purdue  University  Multi-Stage  Axial  Flow 
Research  Compressor. 

Figure  1  presents  a  schematic  of  the  overall  facility,  with 
the  compressor  test  section  shown  in  greater  detail  in  Figures  2 
and  3.  The  three  identical  rotor-stator  stages  are  in  an  annulus 
with  constant  hub  (0.300  m)  and  tip  (0.420  m)  diameters.  All  of 
the  vane  rows  have  41  identical  stators  which  are  mounted  on 
separate  interchangeable  ring  assemblies  which  can  be  indepen¬ 
dently  moved.  The  overall  airfoil  and  compressor  characteristics 
are  presented  in  Table  1. 
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Figure  3. 


Close-Up  View  of  the  Blading  of  the 
Three-Stage  Research  Compressor. 
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TABLE  I 


AIRFOIL  MEAN  SECTION  CHARACTERISTICS  AND 
COMPRESSOR  DESIGN  POINT  CONDITIONS 


Rotor 

Stator 

Type  of  Airfoil 

C4 

C4 

Number 

43 

41 

Chord,  C(nn) 

30 

30 

Solidity,  C/S 

1.09 

1.14 

Aspect  Ratio 

2.0 

2.0 

Thickness/Chord  (X) 

10 

10 

Flow  Rate  (kg/second) 

2.66 

Design  Axial  Velocity  (meters/sec) 

32.0 

Rotational  Speed  (RPM) 

3000 

Number  of  States  3 

Stage  Pressure  Ratio  1.003 

Inlet  Tip  Diameter  (mm)  420 

Hub/Tip  Radius  Ratio  0.715 

Stage  Efficiency,  {%) 


85 
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IV.  PROGRAM  STATUS  AND  RESULTS 


All  of  th#  previously  noted  specific  objectives  established 
for  the  first  year  of  this  research  program  have  been  success¬ 
fully  met. 

*  Dynamic  Instrumentation  of  First  Stage  Vanes 

The  dynamic  instrumentation  for  this  series  of  experiments 
has  been  specified  and  designed.  Por  instrumentation  purposes,  a 
number  of  stator  vanes  have  been  fabricated  from  cast  aluminum. 
In  conjunction  with  Kulite  Corporation,  these  aluminum  airfoils 
have  been  instrumented  along  the  hub,  mean,  and  tip  streamlines 
with  Kulite  thin-line  design  dynamic  pressure  transducers. 

A  structural  vibrations  analysis  has  beet;  performed  on  the 
stator  vanes.  Figure  4  depicts  the  finite  element  representation 
of  the  stator.  The  first  and  second  bend  and  first  torsion  mode 
shapes  and  frequencies  are  presented  in  Figures  5  through  7. 
Based  on  these  results,  a  Campbell  diagram  for  the  first  stage 
vane  row  was  constructed,  Figure  8.  The  43  E  line  represents  the 
first  stage  rotor  blade  excitation  source. 


FREQUENCY:  576.4  Hr.  |SX  BENDING 

Figure  5.  Predicted  Stator  Vane  First  Bend  Mode  Shape 
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FREQUENCY:  3688.4  Hz,  22*  BENDING 

Figure  6.  Predicted  Stator  Vane  Second  Bend  Mode  Shape, 


FREQUENCY:  2607.3  Hz,  |Sl  TORSION 


Figure  7.  Predicted  Stator  Vane  First  Torsion  Mode  Shane 
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*  Dynamic  Data  Acquisition  and  Analysis  System 


I 

! 

In  this  investigation,  both  steady  and  time-variant  data 
will  be  acquired.  Figure  9  presents  a  schematic  of  the  complete 
data  acquisition  and  analysis  system  including  the  interactions 
between  the  steady  and  the  time-variant  systems. 

The  steady-state  data  acquisition  and  analysis  system  is 
centered  around  the  HP-85  computer  seen  in  Figure  10.  These 
steady-state  data  will  define  the  points  of  compressor  operation, 
in  terms  of  overall  and  stage  pressure  ratios  and  corrected  mass 
flow  rate,  at  which  the  unsteady  velocity  and  surface  pressure 
measurements  will  be  obtained. 

The  dynamic  data  acquisition  and  analysis  system  has  been 
specified  and  is  currently  being  made  operational.  This  system 
is  centered  around  a  PDP  11/23  computer  and  an  A-D  multiplexer 
system  capable  of  digitizing  signals  simultaneously  at  rates  of 
up  to  5  megahertz  per  channel,  also  seen  in  Figure  10. 

The  time-variant  quantities  of  fundamental  interest  include 
the  fluctuating  aerodynamic  forcing  function  and  the  resulting 
chordwise  and  spanwise  distributions  of  the  complex  time-variant 
pressure  distributions  of  the  downstream  blade  and  vane  rows. 

The  airfoil  surface  dynamic  pressure  measurements  are  being 
accomplished  with  miniature  high-response  Kulite  thin-line  design 
dynamic  pressure  transducers.  Initially  pairs  of  first  stage 
stator  vanes  are  being  instrumented  such  that  both  the  chordwise 


Figure  9.  Schematic  of  Complete  Data  Acquisition  and  Analysis  System 


Figure  10.  Facility  Computer  and  A-D  Conversion  Systems  . 
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and  spaunwise  tine-variant  pressure  distributions  cam  be  quamtita- 
tively  deternined  in  one  flow  passage.  It  sould  be  noted  that 
the  multi-stage  investigation  will  be  accomplished  by  interchang- 
ing  the  aft  stage  vane  rows  with  the  instrumented  first  stage 
vane  row.  Hence,  it  is  only  necessary  to  instrument  this  first 
stage  vane  row. 

The  t ime-var iant  waUce  measurements  will  be  obtained  by  means 
of  a  calibrated  and  linearized  cross-wire  probe.  The  probe  will 
be  located  at  mid-stator  circumferential  spacing  with  axial  loca¬ 
tion  corresponding  to  mid  rotor -stator  aucial  gap  in  a  passage 
adjacent  to  the  pressure  instrumented  one.  The  probe  will  be 
traversed  randiaLlly,  thereby  defining  the  time-variant  waUce  at  the 
locations  corresponding  to  the  vane  or  blade  surface  spamwise 
dynamic  pressure  measurements.  The  mearn  absolute  flow  angle  will 
be  determined  by  rotating  the  probe  until  a  zero  voltage  differ¬ 
ence  is  obtained  between  the  two  hot-wire  channnels.  This  mean 
arngle  will  then  be  used  as  a  reference  for  calculating  the 
instamtameous  absolute  and  relative  flow  amgles.  The  output  from 
each  channel  will  be  corrected  for  tamgential  cooling  effects  and 
the  individual  fluctuating  velocity  components  parallel  and  nor¬ 
mal  to  the  mean  flow  amgle  calculated  from  the  corrected  quamti- 


ties. 
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*  Data  Analysis  and  Calibration  Techniques 


The  time-variant  data  acquisition  and  analysis  technique 
which  will  be  used  is  based  on  a  data  averaging  or  signal 
enhancement  concept.  The  key  to  such  a  technique  is  the  ability 
to  sample  data  at  a  pre-set  time.  In  this  investigation  the  data 
of  interest  are  being  generated  at  the  blade  passing  frequency. 
Hence,  the  logical  choice  for  the  time  or  data  initiation  refer¬ 
ence  is  the  rotor  shaft.  Por  this  purpose,  am  optical  encoder 
delivering  a  square  wave  voltage  signal  having  a  duration  in  the 
microsecond  range  will  be  used. 

At  each  steady  operating  point,  an  averaged  time-variant 
data  set,  consisting  of  the  hot-wire  and  Kulite  dynamic  pressure 
transducer  signals  averaged  over  V  blade  passages  and  M  rotor 
revolutions,  will  be  obtained.  Each  of  these  digitized  signals 
will  then  be  Fourier  decomposed  into  its  harmonics.  Prom  this, 
both  the  magnitude  of  the  dynamic  data  and  its  phase  lag  as 
referenced  to  the  optical  encoder  pulse  will  be  determined.  To 
then  relate  the  wake  generated  velocity  profiles  with  the  surface 
dynamic  pressures  on  the  instrumented  vanes,  the  rotor  exit  velo¬ 
city  triangles  will  be  examined. 

Figure  11  shows  the  change  in  the  rotor  relative  exit  velo¬ 
city  which  occurs  as  a  result  of  the  presence  of  the  blade.  A 
deficit  in  the  velocity  in  this  relative  frame  creates  a  change 
in  the  absolute  velocity,  as  indicated.  This  velocity  change  is 
measured  with  the  cross  hot  wires.  Prom  this  instantaneous  abso- 


lute  angle  and  velocity,  the  rotor  exit  relative  angle  and  velo¬ 
city  and  the  magnitude  and  phaee  of  the  perturbation  quantities 
will  be  determined. 

As  previously  noted  the  cross  hot-wire  probe  will  be  posi¬ 
tioned  at  mid  rotor-stator  axial  spacing.  To  relate  the  time- 
based  events  as  measured  by  this  hot-wire  probe  to  the  unsteady 
pressures  on  the  vane  surfaces,  the  assumptions  are  made  that: 
(1)  the  wakes  are  identical  at  the  hot-wire  and  the  stator  lead¬ 
ing  edge  planes;  (2)  the  wakes  are  fixed  in  the  relative  frame. 
Figure  12  presents  a  schematic  of  the  rotor  wakes,  the  instru¬ 
mented  vanes,  and  the  hot-wire  probe.  The  rotor  blade  spacing, 
the  vane  spacing,  the  length  of  the  probe,  and  the  axial  spacing 
between  the  vane  leading  edge  plane  and  the  probe  holder  center- 
line  axe  known  quantities.  The  hot-wire  data  will  be  analyzed  to 
yield  the  absolute  flow  angle  and  the  rotor  exit  relative  flow 
angle.  Using  the  two  assumptions  noted,  the  wake  can  be  located 
relative  to  the  hot-wires  and  the  leading  edges  of  the  instru¬ 
mented  van  suction  and  pressure  surfaces.  From  this,  the  times 
at  which  the  wake  is  present  at  various  locations  can  be  deter¬ 
mined.  The  incremented  times  between  occurrences  at  the  hot-wire 
and  the  vane  leading  edge  plane  can  then  be  related  to  phase 
differences  between  the  perturbation  velocities  and  the  airfoil 
surface  pressures. 

The  chordwise  and  spanwise  distributions  of  the  pressure 
difference  on  an  equivalent  single  airfoil  can  then  be  calcu¬ 
lated.  These  data,  as  well  as  the  individual  surface  pressure 


velocity 


Figure  12.  Schematic  of  flow  field  utilized  in  dynamic  data  analysis. 


data,  will  be  normalized  with  respect  to  the  quantity  p'V2*^; 
where  p  is  the  density,  V  is  absolute  velocity,  and  v  the 
transverse  perturbation  velocity  at  the  airfoil  leading  edge. 
These  unsteady  pressure  differential  data  will  then  be  correlated 
with  predictions  obtained  from  an  appropriate  state-of-the-art 
unsteady  aerodynamic  cascade  analysis. 

To  obtain  accurate  time-variant  pressure  data  from  surface 
mounted  Kulite  transducers,  it  is  imperative  that  the  mounted 
transducers  be  dynamically  calibrated  over  the  frequency  range  of 
interest.  Also,  it  is  generally  necessary  to  experimentally 
determine  the  distortion  of  the  unsteady  pressure  signal  associ¬ 
ated  with  the  specific  airfoil  surface  transducer  mounting  confi¬ 
guration.  For  example,  to  minimize  the  disturbance  to  the  air¬ 
foil  contour,  the  transducer  may  be  mounted  from  the  backside  of 
the  airfoil,  sensing  the  unsteady  pressure  via  a  small  volume 
inside  the  airfoil  and  a  static  pressure  tap  on  the  airfoil  sur¬ 
face  of  interest,  as  per  the  airfoil  instrumentation  of  this  pro¬ 
gram.  Another  example  would  include  placing  a  thin  layer  of  RTV 
over  the  flush  mounted  transducer  for  protection.  In  any  event, 
the  distortion  of  the  unsteady  pressure  signal  must  be  experimen¬ 
tally  quantified. 

To  experimentally  determine  the  transfer  function  for  a 
dynamic  pressure  transducer  system,  a  time-variant  pressure  cali¬ 
bration  system  is  required.  Conceptually,  the  simplest  and  most 
direct  calibration  system  involves  an  harmonic  pressure  genera¬ 
tor.  The  calibration  system  of  interest  herein  is  based  on  a 


fluidic  oscillator ,  which  offers  the  potential  of  generating  har¬ 
monic  presssure  signals  at  relatively  high  frequency  values. 

A  fluidic  oscillator,  shown  schematically  in  Figure  13,  is  a 
bi-stable  device  that  results  in  a  fluid  jet  oscillating  between 
two  output  ports.  The  time  duration  of  the  jet  in  each  output 
port,  and  thus  the  oscillator  output  frequency,  is  a  function  of 
both  the  fluid  density  and  the  length  of  the  feedback  line  which 
connects  the  output  and  control  ports.  As  the  fluid  density  and 
the  feedback  line  length  decreases,  the  frequency  of  oscillation 
increases.  It  should  be  noted  that  as  the  frequency  is 
increased,  the  magnitude  of  the  output  signal  is  decreased.  This 
is  not  a  significant  problem  and,  in  fact,  can  be  eliminated  by 
placing  a  fluidic  amplifier  in  series  with  the  fluidic  oscilla¬ 
tor  . 

A  calibration  rig  has  been  designed,  fabricated,  and  util¬ 
ized  to  determine  the  transfer  function  for  two  dynamic  total 
pressure  probe-transducer  systems  with  90°  bends  and  remote 
transducers,  Figure  14,  over  a  frequency  range  up  to  3,000  Hz  and 
correlated  with  predictions  from  the  Bergh-Tij deman  analysis. 
Based  on  the  solution  of  the  linearized,  compressible  Navier- 
Stokes  equations,  this  analysis  predicts  the  dynamic  distortion 
associated  with  a  pressure  tap-tubing-transducer  system. 

Figure  15  presents  the  correlation  of  the  experimentally 
determined  transfer  function  and  the  corresponding  predictions 
for  a  probe-transducer  system  which  has  a  single  resonance  value 


Figure  13.  Schematic  of  Fluidic  Oscillator. 
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in  the  frequency  range  considered.  As  seen,  this  data-analysis 
correlation  is  very  good  when  the  very  small  90°  bend  volume  and 
associated  bend  length  is  considered  in  the  analysis.  However, 
when  this  relatively  small  90°  bend  volume  and  associated  length 
is  neglected,  the  experiment-analysis  correlation  is  less  than 
satisfactory.  Thus,  this  figure  demonstrates  the  viability  of 
the  fluidic  oscillator  based  dynamic  pressure  calibration  system, 
the  validity  of  the  analysis,  and  the  sensitivity  of  the  analysis 
to  relatively  small  variations  in  the  system  geometry. 


*  Energy  Balance  Applied  to  Coupled  Aeroelastic  Response 

The  mathematical  development  of  a  coupled  mode  blade  vibra¬ 
tion  structural  dynamics  model  based  on  an  energy  balance  tech¬ 
nique  has  been  completed.  In  this  technique,  the  energy  input  to 
the  airfoil  system  per  cycle  of  oscillation  is  generated  by  gust 
forces  and  moments  and,  under  certain  conditions,  the  self- 
induced  aerodynamic  forces  and  moments.  The  energy  dissipation 
per  cycle  is  associated  with  the  structural  damping,  the  static 
moment  term  for  coupled  motions,  and  under  certain  conditions, 
the  self-induced  aerodynamic  forces  and  moments. 

The  effects  of  the  various  aerodynamic  parameters  on  the 
coupled  mode  forced  response  of  a  representative  airfoil  have 
also  been  considered.  In  particular,  the  effects  of  the  inlet 
Mach  number,  the  interblade  phase  angle,  the  level  of  structural 
damping,  and  the  cascade  geometry  on  the  coupled  bend ing-tors ion 


aerodynamically  forced  response  of  a  flat  plate  airfoil  cascade 
have  been  investigated. 

A  complete  and  detailed  description  of  this  unique  model  has 
been  presented  as  AIAA  Paper  No.  83-0844  and  published  in  the 
Conference  Proceedings  of  the  24th  AIAA/ASME/ASCE/AHS  Structures, 
Structural  Dynamics  and  Materials  Conference.  This  paper  is 
reproduced  herein  in  the  Appendix. 
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THE  COUPLED  AEROELASTIC  RESPONSE  OF  TURBQMACIIINERY  BLADING 
TO  AERODYNAMIC  EXCITATIONS 

Daniel  lloyniak*  and  Sanford  Fleeter** 

Thermal  Sciences  and  Propulsion  Center 
School  of  Mechanical  Engineering 
Purdue  University 
West  Lafayette,  Indiana  47907 


Abstract 

An  energy  balance  technique  is  developed  to 
predict  the  coupled  bending- tors Ion  mode  aero¬ 
dynamical  ly  forced  response  of  an  airfoil.  In 
this  technique,  the  energy  Input  to  the  airfoil 
system  per  cycle  of  oscillation  is  generated  by 
gust  forces  and  moments  and,  under  certain  condi¬ 
tions,  the  self-induced  aerodynamic  forces  and 
moments.  The  energy  dissipation  per  cycle  is 
associated  with  the  structural  damping,  the  static 
moment  term  for  coupled  motions,  and  under  certain 
conditions,  the  self-induced  aerodynamic  forces 
and  moments.  The  effects  of  the  various  aerodynam¬ 
ic  parameters  on  the  coupled  forced  response  are 
then  considered.  In  particular,  the  effects  of 
the  inlet  Mach  number,  the  interblade  phase  angle, 
the  level  of  structural  damping,  and  the  cascade 
geometry  on  the  coupled  bending-torsion  aerodynam¬ 
ical  1y  forced  response  of  a  flat  plate  airfoil 
cascade  are  demonstrated. 
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Subscripts 

G  gust 

SI  self-induced 

h  translation 

a  torsion 


Superscripts 

I  Imaginary  part 

ft  real  part 

•  derivative  with  respect  to  time 


Introduction 

Aerodynamical ly  Induced  vibrations  of  rotor 
and  stator  airfoils  are  one  of  the  more  common 
sources  of  high  cycle  fatigue  failure  In  gas 
turbine  engines.  Destructive  aerodynamic  forced 
responses  of  fan,  compressor,  and  turbine  blading 
have  been  generated  by  a  wide  variety  of  sources 
Including  igstrean  blades  and/or  vanes,  distortion, 
rotating  stall,  downstream  blades  and/or  vanes, 
surge,  bleeds,  and  random  or  otherwise  unidenti¬ 
fied  sources. 

Failure  level  vibratory  responses  occur  when 
a  periodic  aerodynamic  forcing  function,  with 
frequency  equal  to  a  natural  blade  resonant  fre¬ 
quency,  acts  on  a  blade  row.  The  rotor  speeds  at 
which  these  forced  responses  occur  are  predicted 
with  Campbell  diagrams  which  display  the  natural  . 
frequency  of  each  blade  node  versus  rotor  speed. 
Whenever  these  curves  cross,  aerodynamical ly  In¬ 
duced  forced  responses  are  possible.  However,  no 
accurate  prediction  for  the  amplitude  of  the  re¬ 
sulting  stress  can  currently  be  made. 


analyses.  A  linearized  small  perturbation  gust 
analysis  Is  used  to  predict  the  time-variant  aero¬ 
dynamics  of  the  fixed  non-responding  airfoil  to 
each  harmonic  of  the  disturbance.  A  self-induced 
unsteady  aerodynamic  analysis  wherein  the  airfoils 
are  assumed  to  be  harmonically  oscillating  is  then 
used  to  predict  the  additional  aerodynamic  effect 
due  to  the  airfoil  response.  Reviews  of  the 
state-of-the-art  of  unsteady  aerodynamics  as 
applied  to  turbomachines,  including  gust  and  self- 
induced  unsteady  aerodynamic  analyses,  are  pre¬ 
sented  in  references  3  and  4. 

The  classical  approach  for  the  prediction  of 
the  aerodynamical ly  forced  response  of  an  airfoil 
is  based  on  Newton's  second  law.  The  gust  and 
self- Induced  unsteady  aerodynamic  analyses  are  used 
to  describe  the  harmonic  forces  and  moments  acting 
on  the  airfoil  in  conjunction  with  a  lumped  para¬ 
meter  description  of  the  airfoil  structural  ard 
inertial  properties.  5 

An  alternative  approach  based  on  an  enerqy 
balance  technique  has  been  developed  by  Hoyniak 
and  Fleeter6  to  predict  the  uncoupled  single- 
degree-of- freedom  forced  response  of  an  airfoil. 

In  this  approach,  an  energy  balance  is  conducted 
between  the  unsteady  aerodynamic  work  and  the 
energy  dissipated  by  the  airfoil. 

The  objectives  of  the  research  described  in 
this  paper  include  the  extension  of  the  energy 
balance  technique  to  include  the  more  realistic 
case  of  coupled  bending- tors  Ion  mode  forced  re¬ 
sponse  of  an  airfoil,  and  the  subsequent  demon¬ 
stration  of  the  effects  of  the  various  aerodynamic 
and  structural  parameters  on  this  coupled  response. 


Because  It  is  rarely.  If  ever,  possible  to 
eliminate  all  vibration  excitation  from  the  op¬ 
erating  range  of  turbomachine  blade  raws,  and  not 
possible  to  accurately  predict  the  amplitude  of 
the  stresses  with  present  day  technology,  the 
resonant  stresses  are  not  known  until  the  first 
build-up  test  of  engine  components.  Hhen  exces¬ 
sive  stresses  or  failures  are  then  discovered,  an 
extensive  and  costly  re-design,  re-build,  and 
re-test  program  Is  necessary  In  order  to  effect  a 
viable  solution. 


General  Unsteady  Aerodynamic  Coefficients 

Figure  1  presents  a  schematic  representation 
of  a  two-dimensional  airfoil  section  displaced  in 
both  torsion  and  translation  In  a  uniform  flow  with 
a  superimposed  convected  transverse  sinusoidal 
gust.  Equations  1  and  2  present  the  complex, 
time-dependent  unsteady  lift  and  moment  per  unit 
span,  written  In  influence  coefficient  form,  for 
the  gust  response  and  the  self- Induced  aerodynamic 
cases  respectively. 


The  accurate  first  principles  prediction  of 
the  aerodynamical ly  forced  response  of  a  blade  or 
vane  Involves  the  following  elements.  Spatially 
periodic  variations  In  pressure,  velocity,  and 
flow  direction  of  the  exit  flow  field  of  an  up¬ 
stream  element  appear  as  temporally  varying  In  a 
coordinate  system  fixed  to  the  downstream  row. 

As  a  result.  Individual  airfoils  are  subject  to  a 
time-varying  aerodynamic  forcing  function  which 
can  induce  high  vibratory  stress  levels. 

The  prediction  of  the  aerodynamical ly  forced 
response  vibratory  behavior  of  a  blade  or  vane 
raw  requires  a  definition  of  the  unsteady  forcing 
function  In  terns  of  Its  harmonics.  The  time- 
variant  aerodynamic  response  of  the  airfoil  to 
each  harmonic  of  this  forcing  function  is  then 
assumed  to  be  comprised  of  two  parts.  *  One  Is 
due  to  the  disturbance  being  swept  past  the 
non-responding  fixed  airfoils.  The  second  arises 
when  the  airfoils  respond  to  this  disturbance. 
These  effects  are  modeled  by  means  of  two  distinct 
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H  ■  Hge1"*  describes  the  gust;  (Ah,  Bh), 

(A  ,  B  ),  (A.,  Bj  denote  the  generalized  unsteady 

(x  <1  li  Q 

lift  and  moment  coefficients  due  to  airfoil  trans¬ 
lation,  airfoil  torsion,  and  the  convected  sinu¬ 
soidal  gust,  respectively;  u  Is  the  airfoil  natural 
frequency,  and  w  is  the  gust  forcing  function 
frequency. 
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The  total  unsteady  lift  and  ament  on  the 
airfoil  are  obtained  by  superimposing  the  oust 
response  and  the  self-induced  aerodynamic  forces 
and  aments. 
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Bh  h  ♦  B  «  ]  e’“l 
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It  should  be  noted  that  the  airfoil  response 
occurs  at  the  frequency  of  the  forcing  function. 
Hence  the  gust  frequency*  •»,  has  been  utilized  In 
the  specification  of  the  unsteady  aerodynamic 
forces  and  aments.  Equation  3. 


The  equations  of  notion  describing  the 
coupled  translational  and  torsional  displacement 
of  the  flat  plate  airfoil  depicted  In  Figure  1  are 
given  in  Equation  4. 


*  ♦  v  *4  (1  ♦  1gh)h  «  L(t) 


S  h  ♦  !  «  ♦  I  </  (1  ♦  1g  )«  -  M(t) 
a  a  a  a  a 

where  n  denotes  the  mss  per  unit  span  of  the  air¬ 
foil,  S  Is  the  airfoil  static  noaent  per  unit 
span  abftit  the  elastic  axis,  and  <*a  are  the 
translational  and  torsional  mode  airfoil  natural 
frequencies  of  the  corresponding  undamped  slngle- 
degree-of- freedom  system  respectively,  and  I  is 
the  mss  moment  of  Inertia  per  unit  span  about  the 
elastic  axis.  As  seen,  the  airfoil  structural 
dynamic  system  described  by  Equation  4  Is  strongly 
coupled,  with  this  coupling  associated  with  both 
the  structural  and  the  aerodynamic  characteristics 
of  the  system.  The  structural  coupling  becomes 
significant  whenever  the  aerodynamic  center  does 
not  coincide  with  the  elastic  axis,  thereby  re¬ 
sulting  In  a  non-zero  value  for  the  airfoil  static 
moment.  The  aerodynamic  coupling  arises  because 
the  self-induced  aerodynamic  forces  acting  on  the 
airfoil  are  a  function  of  the  translational  and 
torsional  motions. 

In  this  Investigation,  the  coupled  system  re¬ 
sponse  Is  determined  utilizing  an  energy  balance 
technique.  The  energy  Input  to  the  system  per 
cycle  of  airfoil  oscillation  Is  generated  by  the 
gust  and,  under  certain  conditions,  the  self- 
induced  aerodynamic  forces  and  moments.  6  The 
energy  dissipation  of  the  system  per  cycle  Is 
associated  with  the  system  structural  damping, 
under  certain  conditions  the  self-induced  aero¬ 
dynamic  forces  and  moments,  and  the  static  moment 
term,  $  .  It  should  be  noted  that  for  the  un- 
coupledas1ng1e-degree-of- freedom  case,  the  dis¬ 
sipation  term,  S  ,  Is  not  considered. 
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This  energy  balance  for  the  coupled  transla¬ 
tion-torsion  airfoil  motion  can  be  expressed  as 
follows. 

(WORK)*  ♦  (WORK)]]  ♦  (WORK)!!  «  (W0RK)!) 
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where : 
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(WORK) j 


(WORK)jj,  (WORK)® 


work  done  by  the  self- 
induced  aerodynamic 
forces  and  moments  in 
translation  and  torsion, 
respectively. 

work  done  by  the  airfoil 
translation  and  torsion 
structural  damping, 
respectively. 

work  associated  with  the 
static  moment,  S  ,  in 
translation  and  torsion, 
respectively. 

work  done  by  the  gust 
aerodynamic  forces  and 
moments  in  translation 
and  torsion,  respect¬ 
ively. 


It  should  be  noted  that  the  work  associated 
with  the  self-induced  aerodynamic  forces  and 
moments,  (WORK)jlj  and  (W0RK)|p  have  been  written 
as  dissipation  terms.  However,  under  certain  con¬ 
ditions  they  may  actually  represent  energy  input 
terms.  Mathematically,  this  Is  accounted  for  based 
on  the  sign  Of  these  work  terms.  Also,  matrix 
techniques  will  be  utilized  to  solve  Equation  5  in 
order  to  determine  the  motion  of  the  airfoil 
system. 

The  unsteady  work  done  by  the  self-induced 
aerodynamic  terms,  (WORK)^  and  (W0RK)®j,  over  one 

cycle  of  vibration  can  be  calculated  from  Equation 

6, 


(WORK)*,  «-  j  Re[  L$I(t)  dh] 
(WORK)jj  «  ^  Re[  Mjjtt)  da] 


where  j ( t)  and  Mj,(t)  denote  the  unsteady  lift 

and  moment  associated  with  the  self-induced  har¬ 
monic  motion  of  the  airfoil,  per  Equation  2.  Also, 
the  negative  sign  associated  with  the  aerodynamic 
lift  integral  of  Equation  6  is  a  result  of  the 
standard  aeroelastlc  convention  that  positive  lift 
and  positive  translational  displacement  are  defined 
in  opposite  directions. 

Carrying  out  the  integrations  specified  by 
Equations  2  and  6  results  In  the  following  ex¬ 
pressions  for  the  unsteady  work  done  by  the  self- 
induced  aerodynamic  forces  and  moments. 
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The  angle  xh  denotes  the  phase  angle  between  the 
incasing  gust  and  the  resulting  translational 
notion  of  the  airfoil.  Similarly,  the  angle  X 
is  the  phase  angle  between  the  Incoming  gust  and 
the  airfoil  torsional  motion. 

The  energy  Input  to  the  airfoil  is  determined 
by  calculating  the  work  per  cycle  of  oscillation 
of  the  unsteady  lift  and  moment  associated  with 
the  Incoming  gust. 


(W0RK)|J  -  -  j  Re[L6(t)dh] 
(WORK)q  -  j  Re[KG(t)da] 


Again,  the  negative  sign  is  Included  in  the  lift 
integral  as  a  result  of  the  differences  in  the 
definitions  of  positive  lift  and  positive  trans¬ 
lational  displacement. 

The  evaluation  of  these  Integrals,  Equations 
8  and  2,  yields  the  following  for  the  gust  work. 

(W0RK)g  •  -vpb^hJfAg^COSUj,)  - 
(AG)RSIM(Xh)] 
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The  energy  dissipated  by  the  structural  damp¬ 
ing  over  one  cycle  of  oscillation  can  be  calcu¬ 
lated  as  follows. 


(MOWC)^  »  j  Re[1ghmwjhdh] 
(M0RK)“  -  If  Re[1gaIau,2Qd«3 


Carrying  out  the  Integrations  yields  the 
following: 


The  unsteady  work  associated  with  the  static 
moment,  S^,  merits  some  additional  comment.  This 
is  the  term  which  couples  the  system  equations  of 
motion  when  the  elastic  axis  and  the  aerodynamic 
center  do  not  coincide.  Physically,  the  S  term 
can  be  considered  as  either  a  torsional  force 
applied  to  the  airfoil  as  the  result  of  a  unit 
translation  airfoil  displacement  or  as  a  transla¬ 
tional  force  resulting  from  the  application  of  a 
unit  torsional  displacement.  Thus,  for  an  airfoil 
undergoing  a  translation  mode  oscillation,  there 
is  an  apparent  torsional  load  that  can  either  do 
work  on  or  extract  work  from  the  airfoil  system. 

The  phase  relationship  between  the  translational 
and  torsional  displacement  determines  whether  this 
apparent  torsional  load  adds  or  dissipates  energy 
in  the  airfoil  system. 

The  unsteady  work  per  cycle  of  oscillation 
associated  with  the  static  moment  can  be  determined 
per  Equation  12. 
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The  evaluation  of  these  integrals  leads  to  the 
following. 
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As  previously  noted,  the  sign  of  the  work  associ¬ 
ated  with  the  static  moment  is  determined  by  the 
phase  relationship  between  the  translational  and 
torsional  displacements  per  Equation  13. 

The  substitution  of  the  various  work  expres¬ 
sions,  Equations  7,  9,  11,  and  13,  into  the  energy 
balance.  Equation  5,  yields  the  following  matrix 
expression  for  the  coupled  translational  and  tor¬ 
sional  airfoil  displacements. 
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The  pdas*  angles  between  the  translational 
and  torsional  displacements,  xh  and  x^  respect* 

ively.  can  be  determined  in  various  ways.  For 
example,  they  can  be  determined  directly  from  the 
classical  solution  of  the  coupled  system  equations 
of  motion.  Equation  4.  per  reference  S.  For  con¬ 
venience  and  brevity,  these  phase  angle  relation¬ 
ships  are  presented  In  Equation  15. 


In  tanas  of  the  airfoil  structural  properties  and 
the  unsteady  aerodynamic  forces  and  nonants  in 
the  Nomenclature. 


Unsteady  Aerodynamic  Forces  and  Moments 

L^,  N^,  L#  and  N#  denote  the  standard  fora 

for  the  unsteady  aerodynamic  forces  and  moments 
and  represent  the  unsteady  lift  and  moment  In 
translation  and  torsion,  respectively  calculated 
about  the  airfoil  1/4  chord.  These  standard  form 
1/4  chord  coefficients  are  related  to  the  general 
unsteady  aerodynamic  Influence  coefficients 
through  the  following  relationships. 


<’«  * 

#g  ‘  Si 

(16) 

Bh  "  Mh  *  <1/2  +  #>Lh 

B  »  M  -  (1/2  +  a)(L  +  Mj+O/2  +  a)2Lh 
a  a  an  n 

where:  a  is  the  dimensionless  distance  of  the  tor¬ 
sion  axis  measured  from  the  airfoil  mid-chord,  as 
depicted  in  Figure  1. 


Resul ts 


To  demonstrate  the  effects  of  the  various 
parameters  on  the  coupled  forced  response  of  turbot 
machinery  blading,  the  unsteady  aerodynamic  analy¬ 
sis  of  reference  7  will  be  utilized. 

This  analysis  predicts  the  gust  and  self- 
induced  unsteady  aerodynamic  forces  and  moments  for 
a  cascade  of  flat  plate  airfoils  at  zero  incidence 
in  a  compressible  flow  field.  The  structural  pro¬ 
perties  for  this  parametric  study  are  based  on  an 
airfoil  with  a  5.08  cm  chord,  a  four  percent 
thickness-to-chord  ratio,  and  an  aspect  ratio  of  3. 
This  representative  airfoil  has  a  natural  trans¬ 
lation  mode  frequence,  «h,  12X  larger  than  the 

natural  torsion  mode  frequency,  u»a.  Hence 

uju  »  1.12.  Specific  parameters  to  be  varied 
n  a 

Include  the  structural  damping  value,  the  cascade 
solidity,  the  interblade  phase  angle,  the  inlet 
Mach  number,  and  the  elastic  axis  location.  In 
addition,  to  clearly  demonstrate  the  coupling 
effects,  the  airfoil  stiffness  Is  varied  so  as  to 
make  the  torsional  and  translation  mode  natural 
frequencies  nearly  equal  In  value. 

The  coupled  forced  response  characteristics 
of  the  representative  airfoil  determined  utilizing 
this  energy  balance  technique  are  presented  as  a 
function  of  the  nondimenslonal  ratio  of  the  gust 
forcing  frequency  to  the  airfoil  natural  torsion 
mode  frequency,  3/u^.  Hence,  values  of  1.00  and 

1.12  for  this  nondimenslonal  frequency  ratio  cor¬ 
respond  to  the  situations  where  the  forcing 
function  frequency  Is  equal  to  the  natural  torsion 
and  translation  mode  frequencies,  respectively,  of 
the  representative  airfoil.  Results  presented  in¬ 
clude  the  translational  amplitude,  the  torsion  mode 
amplitude,  and  the  phase  angles  between  the  aero¬ 
dynamic  forcing  function  and  the  resulting  airfoil 
as  a  function  of  the  nondimenslonal  frequency  ratio. 

The  effect  of  structural  damping  on  the 
coupled  forced  response  of  this  representative 
airfoil  is  presented  In  Figures  2  through  5. 

Figures  2  and  3  present  the  translational  and 
torsional  responses  of  the  airfoil  with  structural 
damping  as  a  parameter  when  the  forcing  function 
frequency  is  in  the  neighborhood  of  the  natural 
torsion  mode  frequency,  <ua.  Similarly,  Figures 

4  and  5  present  the  translational  and  torsional 
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airfoil  responses,  respectively,  with  structural 
damping  as  the  parameter  when  the  forcing  function 
frequency  Is  near  to  «h. 

From  these  figures  it  can  be  seen  that  In¬ 
creased  structural  damping  results  in  decreased 
amplitude  of  response.  It  should  be  noted  that 
the  maximum  response  amplitudes  do  not  occur  at 
frequency  ratio  values  of  1.00  or  1.12  even  for 
the  special  case  of  zero  structural  damping.  This 
is  a  result  of  the  self-induced  unsteady  aero¬ 
dynamic  forces  which  correspond  in  this  case  to 
aerodynamic  damping.  Thus,  even  when  there  is  no 
structural  damping  in  the  airfoil  sytem,  the 
self-induced  aerodynamics  can  generate  aerodynamic 
damping. 

Also  as  anticipated,  and  seen  in  Figures  2 
through  S.  when  the  gust  forcing  function  fre¬ 
quency  is  near  to  the  natural  torsion  mode  fre¬ 
quency.  the  torsional  response  amplitude  is  an 
order  of  magnitude  greater  than  the  corresponding 
coupled  translation  mode  response,  Figures  2  and 
3.  Analogously,  when  the  forcing  function  fre¬ 
quency  Is  in  the  near  vicinity  of  the  natural 
translation  mode  frequency,  the  translational 
response  is  an  order  of  magnitude  greater  than  the 
coupled  torsion  mode  response.  Figures  4  and  S. 
Thus,  for  the  other  parameters  considered  herein, 
only  the  coupled  forced  torsional  response  will  be 
presented  for  the  forcing  function  frequency  near 
to  the  airfoil  natural  torsional  frequency 

(u/uAtl),  and  only  the  coupled  forced  transla¬ 
tional  response  for  forcing  function  frequencies 
near  the  natural  translational  frequency 

(u/u  All.  12). 
a 

The  effects  of  cascade  solidity  (C/S)  and 
stagger  angle  on  the  coupled  forced  response  amp¬ 
litudes  are  shown  In  Figures  6  and  7  and  Figures 
8  and  9,  respectively.  As  seen,  both  the  coupled 
torsional  and  translational  forced  response  amp¬ 
litudes  Increase  with  Increasing  values  of  the 
solidity  and  the  stagger  angle.  The  variations  In 
the  solidity  result  In  a  somewhat  more  pronounced 
effect  on  the  translational  response  than  on  the 
torsional  response.  Figures  6  and  7.  However,  the 
stagger  angle  variation  has  a  greater  effect  on 
the  torsional  response  than  on  the  translational 
one,  per  Figures  8  and  9.  Also,  variations  In  the 
stagger  angle  have  a  much  greater  effect  on  the 
response  amplitudes,  both  torsional  and  transla¬ 
tional,  than  do  the  variations  in  the  solidity. 

It  should  be  noted  that  the  design  trend  for 
modern  compressors  Includes  Increased  solidity  and 
stagger  angle  values.  Hence,  these  results  Indi¬ 
cate  that  Increased  forced  response  problems,  pri¬ 
marily  associated  with  Increased  stagger  angles, 
■right  be  anticipated. 

Figures  10  and  11  show  the  effect  of  Inter- 
blade  phase  angle  on  the  coupled  airfoil  forced 
response  characteristics.  As  seen  the  effect  of 
variations  In  this  parameter  on  the  coupled  tor¬ 
sional  forced  response  are  quite  different  than 
Its  effect  on  the  coupled  translational  response. 
For  the  torsion  mode.  Figure  10,  there  Is  only  an 
extreme  maximum  response  amplitude  for  a  0°  Inter¬ 
blade  phase  angle  value.  For  all  other  values 
considered,  the  response  amplitude  curves  are  re¬ 
latively  flat  and  decreased  from  the  0°  case  by 
approximately  759  to  859.  The  coupled  transla¬ 


tional  response  amplitudes  all  exhibit  an  extreme 
maximum  response  amplitude,  l.e.  none  of  the  amp¬ 
litude  response  curves  are  relatively  flat.  Figure 
11.  Also,  although  an  interblade  phase  angle  value 
of  0°  results  In  the  maximum  forced  translational 
response  amplitude,  varying  the  interblade  phase 
angle  away  from  0°  results  In  noticably  smaller 
decreases  in  this  maximum  amplitude  than  were  noted 
for  the  forced  torsional  response.  In  particular, 
for  both  the  torsional  and  the  translational  cases, 
the  largest  decrease  in  the  maximum  response  amp¬ 
litude  Is  associated  with  changing  the  interblade 
phase  angle  value  from  0°  to  -90°.  However,  the 
maximum  forced  translational  response  decreased 
only  459  whereas  the  corresponding  forced  torsional 
response  decreased  by  859.  In  addition,  as  the 
smallest  coupled  torsional  and  translational  re¬ 
sponse  amplitudes  correspond  to  negative  interblade 
phase  angle  values.  In  particular  -90°,  then  this 
would  appear  to  be  a  possibly  desirable  forced  re¬ 
sponse  design  condition.  In  terms  of  the  rotor- 
stator  Interaction  forced  response  problem,  a 
negative  Interblade  phase  angle  value  corresponds 
to  a  backward  traveling  wave  as  viewed  from  the 
stator  vane  frame  of  reference  and  arises  when  the 
number  of  rotor  blades  Is  greater  than  the  number 
of  stator  vanes. 

The  effect  of  the  Inlet  flow  Mach  number  on 
the  coupled  torsional  and  translational  forced  re¬ 
sponse  characteristics  of  the  representative  air¬ 
foil  are  demonstrated  In  Figures  12  and  13, 
respectively.  Increasing  the  Mach  number,  which 
also  corresponds  to  decreasing  the  reduced  fre¬ 
quency  For  a  specific  gust  disturbance.  Increases 
both  the  coupled  torsional  and  translational  re¬ 
sponses,  with  a  very  dramatic  Increase  associated 
with  an  increase  In  Mach  number  from  0.6  to  0.8. 
Also,  for  Mach  numbers  of  0.6  and  smaller,  the  in¬ 
crease  in  the  maximum  response  amplitude  is  much 
larger  for  the  torsional  mode  than  for  the  trans¬ 
lational  mode.  In  addition,  the  forcing  function 
frequency  at  which  the  maximum  amplitude  response 
occurs  is  hiqhly  dependent  on  the  Mach  number  for 
the  torsional  response,  but  Is  almost  Independent 
of  this  forcing  frequency  for  the  translational 
response  case. 

The  effect  of  moving  the  elastic  axis  loca¬ 
tion  to  the  1/4-chord  and  the  3/4-chord  locations 
on  the  coupled  forced  response  of  the  representa¬ 
tive  airfoil  is  presented  In  Figures  14  through  17. 
These  results  Indicate  that  not  only  the  response 
amplitude  but  also  the  forcing  function  frequency 
at  which  the  maximum  amplitude  response  occurs  is 
strongly  dependent  on  the  position  of  the  elastic 
axis.  The  maximum  torsional  and  translational  re¬ 
sponse  amplitudes  are  found  with  the  elastic  axis 
at  the  3/4-chord  location.  Shifting  the  elastic 
axis  forward  to  the  1/4-chord  location  results  In 
significantly  decreasing  both  the  torsional  and 
translational  response  amplitudes  near  the  airfoil 
torsional  natural  frequency  but  only  slightly  de¬ 
creasing  these  response  amplitudes  near  the  airfoil 
translational  natural  frequency.  Also,  Figures  14 
and  15  demonstrate  that  shifting  the  elastic  axis 
location  aft  from  mid-chord  to  3/4-chord  results 
in  the  translational  response  In  the  neighborhood 
of  the  natural  torsional  frequency  (u/u^adl.o) 

being  of  the  same  order  of  magnitude  as  the  tor¬ 
sional  response  amplitude.  This  Is  an  Indication 
of  the  significant  additional  coupling  that  arises 
when  the  elastic  axis  and  the  aifoil  center  of 
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gravity  do  not  coincide.  Similar  results  are 
noted  when  the  forcing  function  frequency  Is  near 
to  the  natural  translational  frequency.  Figures 
16  and  17. 

To  demonstrate  the  Interaction  of  two  very 
closely  spaced  torsional  and  translation  modes, 
the  bending  stiffness  of  the  representative  air¬ 
foil  was  altered  such  that  the  ratio  of  the  trans¬ 
lation  to  torsion  natural  frequencies  became  1.02. 
The  variation  with  Interblade  phase  angle  values 
of  these  closely  coupled  forced  torsional  and 
translational  responses  was  then  Investigated  with 
the  results  presented  in  Figures  18  and  19.  The 
coupled  torsional  response  results.  Figure  18, 
show  that  while  this  response  at  a  frequency  ratio 
of  0.990  decreases  as  the  Interblade  phase  angle 
Is  varied  from  0°  to  ■*•90°,  -90°,  and  180s,  this 
response  Increases  by  up  to  65%  at  a  frequency 
ratio  of  1.015. 

Figure  19  presents  the  variation  of  the 
forced  translational  response  with  the  Interblade 
phase  angle  as  the  parameter.  In  contrast  to  the 
torsional  response  results  of  Figure  18,  minimal 
translational  response  Is  noted  when  the  forcing 
function  frequency  Is  near  the  airfoil  natural 
torsional  frequency.  However,  when  the  forcing 
function  frequency  Is  near  the  airfoil  natural 

translational  frequency  (u/u  «  1 .015),  the  forced 

a 

translational  amplitude  Is  significant  and  in¬ 
creases  as  the  Interblade  phase  angle  value  Is 
varied  from  0  to  180°.  A  comparison  of  Figures 
18  and  19  at  a  frequency  ratio  near  1.015  reveals 
a  correspondence  between  the  forced  torsional  re¬ 
sponse  and  the  translational  one.  This  Is  caused 
by  the  high  degree  of  coupling  that  exists  between 
'the  translation  and  torsional  modes  for  this 
altered  airfoil,  even  when  the  elastic  axis  coin¬ 
cides  with  the  airfoil  mid-chord. 


Suwaary  and  Conclusions 

In  this  Investigation,  the  coupled  struc¬ 
tural  dynamic  airfoil  system  forced  response  was 
determined  utilizing  an  energy  balance  technique. 
The  energy  Input  to  the  system  per  cycle  of  oscil¬ 
lation  was  generated  by  the  gust  and,  under  certain 
conditions,  the  self-induced  aerodynamic  forces 
and  moments.  The  energy  dissipation  per  cycle  of 
oscillation  was  associated  with  the  system  struc¬ 
tural  damping,  under  some  conditions  the  self- 
induced  aerodynamic  forces  and  moments,  and  the 

static  moment  term  S  .  It  should  be  noted  that 
(1 

for  the  uncoupled  s 1 ngle- deg ree-of- freedom  case, 
the  dissipative  coupling  mechanism,  Sa,  Is  not 

considered. 

The  effects  of  the  various  aerodynamic  para¬ 
meters  on  the  coupled  translational  and  torsional 
mode  forced  response  of  a  representative  airfoil 
was  then  considered.  This  study  showed  the  In¬ 
creased  coupling  between  the  torsion  and  trans¬ 
lation  modes  as  the  corresponding  undamped  natural 
frequencies  approach  one  another.  It  was  also 
demonstrated  that  the  coupled  torsional  and  trans¬ 
lational  forced  response  amplitudes  of  a  repre¬ 
sentative  airfoil  Increased  with: 


*  increased  stagger  angles 

*  increased  Inlet  flow  Mach  numbers 

*  interblade  phase  angle  values  cor¬ 
responding  to  forward  traveling  waves 
for  the  rotor-stator  interaction  case 

*  shifting  of  the  elastic  axis  location 
aft. 

It  should  be  noted  that  variations  in  the  above 
parameters  did  not  affect  the  magnitude  of  the 
resulting  forced  response  equally  nor  did  they 
always  have  equivalent  effects  on  the  torsional 
and  translational  response  modes.  For  example, 
increasing  the  stagger  angle  and  shifting  the  elas¬ 
tic  axis  aft  resulted  in  significantly  larger  in¬ 
creases  in  all  of  the  response  amplitudes  than  did 
increasing  the  solidity.  Also,  the  Interblade 
phase  angle,  for  example,  had  a  much  greater  effect 
on  the  coupled  torsion  mode  responses  than  on  the 
translational  ones.  In  addition,  the  forcinq 
function  frequencies  at  which  the  maximum  torsional 
and  translational  responses  occurred  did  not  gen¬ 
erally  correspond  to  the  airfoil  natural  torsional 
or  translational  frequencies.  This  was  because 
the  self-induced  aerodynamic  forces  and  moments 
result  in  an  aerodynamic  damping  effect  analogous 
to  structural  damping. 
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Figure  7.  Cascade  Solidity  Effect  on  the  Coupled 
Translational  Response. 
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Figure  9.  Effect  of  Stagger  Angle  on  the  Coupled  Figure  11.  Effect  of  Interblade  Phase  Angle  on 

Translational  Response.  the  Coupled  Translation  Mode  Response. 


Figure  10,  Effect  of  Interblade  Phase  Angle  on  Figure  12.  Inlet  Mach  Humber  Effect  on  the  Coupled 
the  Coupled  Torsion  Mode  Response.  Torsional  Response. 
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Figure  13.  Inlet  Mach  Number  Effect  on  the  Coupled  Figure  15.  Effect  of  Elastic  Axis  Location  on  the 
Translational  Response.  Coupled  Translational  Response  for 
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Figure  14.  Effect  of  Elastic  Axis  Location  on  the 
Coupled  Torsional  Response  for  (Z/u  ) 
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Figure  16.  Effect  of  Elastic  Axis  Location  or.  th 
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Figure  17.  Effect  of  Elastic  Axis  Location  on  the 
Coupled  Translational  Response  for 
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Figure  19.  Interblade  Phase  Angle  Effect  on  the 
Coupled  Translational  Response  for 
Nearly  Equal  Natural  Frequencies. 
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Figure  IS.  Interblade  Phase  Angle  Effect  on  the 
Coupled  Torsional  Response  for  Nearly 
Equal  Natural  Frequencies. 
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